Potential energy topology and relaxation processes in a model glass 
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We use computer simulation to investigate the topology of the potential energy F({R}) and to 
search for doublewell potential's (DWP) in a model glass . By a sequence of Newtonian and dissi- 
pative dynamics we find different minima of F({R}) and the energy profile along the least action 
paths joining them. At variance with previous suggestions, we find that the parameters describ- 
ing the DWP's are correlated among each others. Moreover, the trajectory of the system in the 
3N-d configurational phase space follows a quasi- 1-d manifold. The motion parallel to the path is 
characterized by jumps between minima, and is nearly uncorrelated from the orthogonal, harmonic, 
dynamics. 



PACS Numbers : 61.43. Fs, 63.50.+X, 82.20.Wt 

Among many unanswered questions concerning the dy- 
namics of topologically disordered materials, two are re- 
lated to the low energy excitations. The first one regards 
the excitations responsible for the low temperature ther- 
mal properties of glasses. Since the pioneering works of 
Phillips 0] and Anderson et al. who introduced the 
Two Level Systems (TLS) as the origin of the anoma- 
lous specific heat behaviour below 1 K, many works were 
devoted to identifying the degrees of freedom associated 
to the TLS in real glasses, but a definite answer has not 
yet been reached. The second question regards the exci- 
tations responsible for: (z) the thermal behavior in the 
1-^10 K region, where a plateau in the thermal conduc- 
tivity and a peak in c p /T 3 are observed ||; and (ii) the 
boson peak appearing in the inelastic neutron scatter- 
ing (INS) (§-§,[ll) and Raman spectra [0-0. This sec- 
ond question is also largely unanswered, and different 
hypotheses on the nature of these excitations have been 
proposed, ranging from an excess of acoustic-like modes 
to highly anharmonic and localized excitations. A link 
between the above questions is suggested by the so called 
Soft Potential Model (SPM) O, which assumes the ex- 
istence of few higly anharmonic degrees of freedom with 
a potential energy function described by a fourth order 
polynomial. The coefficients of the polynomial are sup- 
posed to be uncorrelated, and to give rise to a large va- 
riety of curves, including doublewell potentials (DWP's). 
Among the DWP's, those that are isolated, have low 
asymmetry and have barrier height resulting in a tun- 
neling splitting around 1 K, are candidate to be TLS. 

Computer Simulations (CS) of model glasses might 
help to answer the previous questions, and, in partic- 
ular, to verify the existence of TLS and to clarify the 
nature of DWP in glasses. To our knowledge the di- 
rect CS inspection of the passage from minimum to min- 
imum in glasses is limited to a few cases. Among them 
we mention the search for TLS performed by Stillinger 
and Weber @ and by Heuer and Silbey The 
conclusions of these works can be summarized as fol- 
lows: i) low barrier bistable degrees of freedom exist in 



binary glasses with considerable variation in their phys- 
ical parameters (barrier height, asimmetry, minima dis- 
tance, etc.); ii) these parameters are strongly correlated, 
DWP with large asymmetry typically have large barrier 
heights; Hi) in agreement with SPM, the coefficients of 
the fourth order polynomial describing the potential en- 
ergy of the reaction coordinate, are found to be uncorre- 
lated; iv) using a factorized probability distribution func- 
tion for these coefficients, it is possible to try an estimate 
of the number of TLS per atom; v) typically from 5 to 
10 atoms participate in the motion in the DPW. 

In this letter we investigate by CS the topology of the 
potential energy hypersurface of a Lennard-Jones glass 
with the aim to check whether the TLS can be observed 
in this system, the characteristics of the DWP's, and the 
extent to which the harmonicity of the system is affected 
by the presence of DWP's. We will discuss the procedures 
used to search the minima of y({R}) and to identify the 
Reaction Coordinate (RC) and the Least Action Path 
(LAP) joining pairs of minima. The main conclusions of 
the present work arc the following: a) no isolated double- 
wells are found in the system, i.e. in the investigated 
temperature range (T > 5 K) all the minima pertain 
to a network, and the thermally activated jumps among 
them are not controlled only by the energetic barriers, 
because, depending on temperature, a significant role is 
also played by the entropic term . b) At variance with 
previous findings [i~^Jl5[| (see item Hi) above) and with 
the hypothesis at the base of the SPM [[l2| the coefficients 
of the polynomial representation of the energy profile felt 
by the RC, are found to be higly correlated, c) The 
degrees of freedom orthogonal to the LAP feel an almost 
harmonic potential, whose curvature is independent on 
the RC, i. .e. on the position along the LAP itself. 

We investigated a sample of N = 864 atoms interact- 
ing via the 6-12 Lennard-Jones potential, here we use 
e/Ks = 125.2 K and a — 0.3405 nm, appropriate for 
Argon. A microcanonical molecular dynamics simulation 
JlTj is carried out at the fixed density of 42 mol/dm 3 , and 
at different temperatures T (T — 15, 12, 10, 7.5 and 6 



1 



K). Each run consists of a sequence of newtonian and dis- 
sipative dynamics. The newtonian trajectory is followed 
for 1 ps (50 integration time steps), subsequently a mod- 
ified steepest descent method procedure (relaxation-like 
dynamics) is applied to quench the system and to find 
the inherent configuration {X,} (i = 1...N). correspond- 
ing to a local minimum. The newtonian dynamics is 
then restarted from the same point in phase space where 
it was interrupted before the quenching, and the proce- 
dure is repeated up to 25,000 times for each temperature. 
The two minimum configurations obtained with succes- 
sive quenching, {X°} and {X*}, are often the same. 

The adopted procedure ensures that different minima, 
when found, are close to each other in the 3N-d config- 
urational space. About 130 different minima have been 
found in the present work: their energies are spread in 
an interval of about 1 J/molc centered at about 300 
J/mole above the absolute minimum representing the 
crystal. All the minima belong to an inter-connected net- 
work and the system is observed to perform temporarily 
closed loop in the configurational phase space. In or- 
der to study the topopogy of the potential energy, and 
to find the easiest way to jump among neighbourING 
minima, a fast and efficient algorithm is developed to 
evaluate the Least Action Path (LAP). The LAP is de- 
termined as the path between X" and X^ that minimize 

the classical action integral, J xa dsy / V(R(s)) — Vq, with 
V Q = mm{V({X?}),V({X£})}. We studied « 100 pairs 
of minima, evaluating the LAP joning them; in Fig. la 
we compare the potential energy profile evaluated along 
the LAP with that along the straight path. Although the 
LAP is not very far from being straight in the configu- 
rational space, as shown in Fig. lb, the potential energy 
along the LAP is significantly lower, making the LAP it- 
self highly preferred for the jump among minima. More- 
over, the present analysis shows that the derivation of the 
DWP properties from the profile of the potential energy 
along the straight path between minima gives unreliable 
results. The potential energy profile along the LAP be- 
tween the minima a and b has been characterized by the 
energy difference between the two minima AV^™, their 
euclidean distance D a b, the displacement of the atom 
that moves most d a b = ma,Xi{\~X.f — X^|}, the barrier 
height AVj^ (measured with respect to the higher mini- 
mum) and finally by the number of participation, defined 
as P ab = (£, |X? - X|| 2 ) 2 / Ei IX? - X?| 4 . It turns out 
that: i) the participation is always around 20; ii) among 
the 20 atoms involved in the jump, 1 or 2 account for 
90% of the entire distance; Hi) for each pair of minima 
AV m is always higher than AV B . 

Next we analyse the statistical correlation among the 
DWP's parameters. In Fig. 2a we see that the dif- 
ferent pairs of parametres that characterized the DWP 
are strongly correlated. To give a statistical signifi- 
cance to the correlation, for each pair of parameters we 



calculated the correlation coefficient r ]T^| . The val- 
ues of r, reported in the first part of Table I, indi- 
cate statistical correlation between the parameters (with 
nDWP = 141, at the 5% of significance level, the thresh- 
old for correlation is r = 0.165). Two other representa- 
tions of the DWP are reported in the literature, these 
are the "u" representation, used in the soft potential 
model @, E = e u [u 2 {x/a) 2 - u 3 (x/a) 3 + (x/a) 4 }, 
and the "w" representation introduced in ref. [ pif , E = 
e [w 2 (x/a) 2 - w 3 (x/cr) 3 +w 4 (x/a) 4 }. The SPM assumes 
that there is no correlation among the parameters of the 
set {u ,u 2 ,u 3 }, while the authors of ref. p5| |, noticing 
the correlation of the u set, made use of the uncorre- 
lation among {w 2 , W3, W4} to determine the number of 
TLS in model glasses. As can be seen in Figs. 2b and 2c, 
and in Table I, we found statistical correlation also in the 
case of the w set. It is beyond the scope of the present 
work to explain this inconsistency, that can lie on the 
different glass examined or, more likely, on the different 
procedure used to determine the potential energy profile 
(LAP here, straight path in |jl5|). 

In order to describe the low temperature anomalies of 
glasses, the most important quantity is the total split- 
ting of the ground state associated to the TLS. To our 
knowledge, only one attempt has been made to estimate 
the energy splitting due to the tunneling through TLS 
energy barriers. It a) assumes the existence of isolated 
pairs of minima and b) treats the problem as it were 
1-dimensional (1-d) From the present simulation, 

hypothesis a) can be neither cibfirmed nor rejected be- 
cause we did not directly identify any candidate TLS 
(i. e. pairs of minima with asymmetry < 1 K); recent re- 
sults on Argon clusters ]l9| ] indicate that pairs of minima 
with low asymmetry are "isolated" in the sense that a 
third quasi-degenerate, adjacent minimum has not been 
observed. As for item b), the (3N — 3)-dimensional tun- 
neling problem can be treated as an 1-dimensional one if 
the Schrodinger equation can be factorized into 3iV — 3 
independent equations p0(| . We shall assume that the 
reduction to 1— D can be made under the less restric- 
tive condition that the relevant classical path (i. e. the 
LAP) is independent of all the others. To this end we 
have studied the curvatures of the potential energy sur- 
face along the LAP by calculating the dynamical matrix 
eigenvalues \j(n) (j = 1, ..,3iV — 3) and the correspond- 
ing eigenvectors in 42 (n = 1..42) equally spaced config- 
urations along the LAP itself. In Fig. 3, we report, for 
each configuration along the LAP joining a typical pair of 
minima, the 20 (j = 1..20) lowest frequencies uoj = | 
(we have assigned the minus sign to the frequencies as- 
sociated to negative eigenvalues) . As can be seen, in the 
present case only one eigenvalue becomes negative, indi- 
cating a first order saddle point in the 3N-dimensional 
space, and only the lowest frequencies change apprecia- 
bly along the LAP. In Fig. 4 we report, for the seven 
lowest eigenvalues, the projections of the eigenvectors on 
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the local tangent at the LAP. Along the path, the main 
contribution comes from the lowest eigenvalue, while ap- 
proaching the minima, an increasing contribution comes 
from other eigenvalues. The same results are found for all 
pairs of minima. They indicate that the dynamics paral- 
lel and orthogonal to the reaction coordinate are nearly 
independent: quasi-harmonic vibrations control the or- 
thogonal dynamics, by a set of eigenfrequencies that are 
independent of the position along the path. The reaction 
coordinate is associated with a single eigenvalue (often 
the lowest one), as shown by the large value of the pro- 
jection of only one eigenvector in the path direction. 

In conclusion, the picture emerging from the depicted 
scenario is that of a network of connected minima, where 
each pair is joined together by a 1-dimensional path. The 
reaction coordinate along the path follows a LAP in the 
3N-configurational space, described by DWP with statis- 
tically correlated parameters (asymmetry, height, etc.). 
The potential energy experienced by the system along 
the LAP is significantly lower than that one along other 
paths, like, for example, the straight path. Orthogonally 
to the RC, the dynamics is harmonic, and the set of fre- 
quencies are nearly independent on the specific value of 
the RC. In the present simulation a single interconnected 
quasi-l-dim path has been found, indicating that no more 
than one DWP is active at the time. Extrapolating this 
result, we can state that the number of DWP 's coexisting 
with harmonic excitations is Ndwp < 10~ 3 atoms . 
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TABLE I. The upper right part of the 
ports the correlation coefficients r among 
parameters used to describe the DWP in 
resentation {AV"™, AV B , D, d, P}, and 
ent polynomial representations {ui2, W3, 
us}. The lower left part of the correlation 
whether the two parameters are (»)or are 
at the 5% level of confidence. 
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CAPTIONS 



Fig. 1 - a) Comparison between the potential energy profile along 
the LAP (•) and along the straight path (o) joining two 
typical minima. E x is the energy of the crystalline minima, 
b) Orthogonal distance between the LAP and the straight 
path as a function of the reaction coordinate. 

Fig. 2 - Example of correlation between the parameters describing 
the DWP. Each point represent the values found for the 
couple of parameters a) (AV m , D), b) (w4, 11)3), and c) (113, 
112) in the description of each DWP. 



Fig. 3 - The 20 lowest eigen-frequencies u> = ^/|A| of the system 

evaluated in 42 cquispaccd atomic configurations along the 
LAP are reported as a function of the reaction coordinate 
for a typical pair of minima. Negative eigenvalues A arc 
reported as negative eigen-frcquencies cj. 

Fig. 4 - The projection on the LAP of the six lowest frequency 
cigen- vectors of the modes of Fig. 3 are reported as a func- 
tion of the reaction coordinate. 
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Fig. 1 - "Potential energy topology "by F. Demichelis et al. 
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Fig. 2 - "Potential energy topology " by F. Demichelis et al. 
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Fig. 4 - "Potential energy topology " by F. Demichelis et al. 



